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Abstract. We introduce an (evolution) algebra identifying the coeffi- 
cients of inheritance of a bisexual population as the structure constants 
of the algebra. The basic properties of the algebra are studied. We 
prove that this algebra is commutative (and hence flexible), not associa- 
tive and not necessarily power associative. We show that the evolution 
algebra of the bisexual population is not a baric algebra, but a dibaric 
algebra and hence its square is baric. Moreover, we show that the al- 
gebra is a Banach algebra. The set of all derivations of the evolution 
algebra is described. We find necessary conditions for a state of the 
population to be a fixed point or a zero point of the evolution operator 
which corresponds to the evolution algebra. We also establish upper es- 
timate of the limit points set for trajectories of the evolution operator. 
Using the necessary conditions we give a detailed analysis of a special 
case of the evolution algebra (bisexual population of which has a pref- 
erence on type "1" of females and males). For such a special case we 
describe the full set of idempotent elements and the full set of absolute 
nilpotent elements. 



1. Introduction 

The action of genes is manifested statistically in sufficiently large com- 
munities of matching individuals (belonging to the same species). These 
communities are called populations [13]. The population exists not only in 
space but also in time, i.e. it has its own life cycle. The basis for this phe- 
nomenon is reproduction by mating. Mating in a population can be free or 
subject to certain restrictions. 

The whole population in space and time comprises discrete generations 
Fq,Fx, . . . The generation F n+ \ is the set of individuals whose parents be- 
long to the F n generation. A state of a population is a distribution of 
probabilities of the different types of organisms in every generation. Type 
partition is called differentiation. The simplest example is sex differentia- 
tion. In bisexual population any kind of differentiation must agree with the 
sex differentiation, i.e. all the organisms of one type must belong to the 
same sex. Thus, it is possible to speak of male and female types. 

The evolution (or dynamics) of a population comprises a determined 
change of state in the next generations as a result of reproduction and selec- 
tion. This evolution of a population can be studied by a dynamical system 
(iterations) of a quadratic stochastic operator. 
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The history of the quadratic stochastic operators can be traced back to the 
work of S. Bernstein pp. During more than 85 years this theory developed 
and many papers were published (see e.g. p^[6ti9 | fT3 ].fT5HT8] ) . In recent years 
it has again become of interest in connection with numerous applications to 
many branches of mathematics, biology and physics. 

A quadratic stochastic operator (QSO), V, has meaning of a free popula- 
tion evolution operator, which arises as follows: Consider a free population 
consisting of m species. Let x° = (a^, . . . , x^) be the probability distri- 
bution of species in the initial generations, and Pij t k the probability that 
individuals in the ith and jth species interbreed to produce an individual 
k. Then the probability distribution x' = (x^, . . . ,x' m ) of the species in the 
first generation can be found by the total probability i.e. 

m 

(1.1) x' k = ^2 p ij,k^iXj, k = 1, . . . , m 

where the cubic matrix P = P(V) = (Pij^)^ k=i satisfies the following 
conditions 

TO 

(1.2) Pij, k > 0, Y. Pi ^ k = X > G i 1 ' • • • ' m >- 

k=l 

This means that the association x° i— )■ x' defines a map V of the simplex 

m 

(1.3) S m ~ l ={x = (xi, . . . , x m ) € M m : x< > 0, x { = 1} 

i=l 

into itself, called the evolution operator. 

The population evolves by starting from an arbitrary state x^ , then pass- 
ing to the state x' = V(x^), then to the state x" = V(V(x^)), and so on. 

There are many papers devoted to study of the evolution of the free 
population, i.e. to study of dynamical system generated by the QSO (jl.ip . 
see e.g. PEESIEHIB] . 

There exist several classes of non-associative algebras (baric, evolution, 
Bernstein, train, stochastic, etc.), whose investigation has provided a num- 
ber of significant contributions to theoretical population genetics. Such 
classes have been defined different times by several authors, and all algebras 
belonging to these classes are generally called genetic. Etherington intro- 
duced the formal language of abstract algebra to the study of the genetics 
in his series of seminal papers [HHS] . In recent years many authors have tried 
to investigate the difficult problem of classification of these algebras. The 
most comprehensive references for the mathematical research done in this 
area are (EHHlIISlESg . 
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In [13] an evolution algebra A associated to the free population is intro- 
duced and using this non-associative algebra many results are obtained in 
explicit form, e.g. the explicit description of stationary quadratic operators, 
and the explicit solutions of a nonlinear evolutionary equation in the absence 
of selection, as well as general theorems on convergence to equilibrium in 
the presence of selection. 

In [19] a new type of evolution algebra is introduced. This algebra also 
describes some evolution laws of genetics and it is an algebra E over a 
field K with a countable natural basis e\,e2,--- and multiplication given 
by e-i&i = Ylj a ij e j> e i e j = if £ ^ j. Therefore, is viewed as "self- 
reproduction" . 

In this paper we consider a bisexual population (BP) and define an evo- 
lution algebra (EA) using inheritance coefficients of the population. This 
algebra is a natural generalization of the algebra A of free population. The 
evolution algebra of a bisexual population (EABP) is different from the EA 
defined in pi]. 

The paper is organized as follows. In Section 2 we give evolution opera- 
tor of BP. Section 3 contains the definition of EABP. Section 4 is devoted 
to basic properties of the algebra and therein we prove that the EABP is 
commutative but not associative and not power associative. In Section 5 we 
show that the EABP is not a baric algebra. In Section 6 we prove that the 
EABP is a dibaric algebra, hence its square is a baric algebra. In Section 7 
we prove that the EABP is a Banach algebra. In Section 8 we describe the 
set of all derivations of EABP. Section 9 is devoted to study dynamics of 
the evolution operator, which corresponds to the evolution algebra. We find 
necessary conditions for a state of the population to be a fixed point or a 
zero point of the evolution operator. We also establish upper estimate of the 
limit points set for trajectories of the evolution operator. In the last section 
we give a detailed analysis of a special case of the evolution algebra (bisexual 
population of which has a preference on type "1" of females and males). For 
such a special case we describe the full set of idempotent elements and the 
full set of absolute nilpotent elements. 



In this section following [13] , we describe the evolution operator of a BP. 
Assuming that the population is bisexual we suppose that the set of females 
can be partitioned into finitely many different types indexed by {1, 2, . . . , n} 
and, similarly, that the male types are indexed by {1, 2, . . . , v}. The number 
n + v is called the dimension of the population. The population is described 
by its state vector (x,y) in S 1 " -1 x S v ~ x , the product of two unit simplexes 
in M n and W respectively. Vectors x and y are the probability distributions 
of the females and males over the possible types: 



2. Evolution operator of a BP 



n 



V 



Xi > 0, 



Xi = 1; yi> 0, = 1. 
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Denote S = S n ~ x S u . We call the partition into types hereditary if 
for each possible state z = (x,y) £ S describing the current generation, the 
state z' = (x',y ! ) € S is uniquely defined describing the next generation. 
This means that the association z i— >■ z 1 defines a map V : S — > S called the 
evolution operator. 

For any point z^ € S the sequence z^ = V(z^ t ~ 1 ^), t = 1,2,... is called 

the trajectory of z(°\ 
(f) (rri) 

Let ■ and P ifc { be inheritance coefficients defined as the probability 
that a female offspring is type j and, respectively, that a male offspring is 
of type I, when the parental pair is ik (i,j = 1, . . . , n; and k, I = 1, . . . , v). 
We have 

n v 
3=1 1=1 

Let z' = (x' , y') be the state of the offspring population at the birth stage. 
This is obtained from inheritance coefficients as 



n,u n,u 

(2-2) 4 = £ pw*^ y ; = £ 

i,fc=l j,fc=l 

We see from (|2.2p that for a BP the evolution operator is a quadratic 
mapping of S into itself. But for free population the operator is quadratic 
mapping of the simplex into itself given by (jl.ip . 



3. Definition of the EABP 

In this section we give an algebra structure on the vector space 
which is closely related to the map (|2.2p . 

Consider {ei, . . . , e n+u } the canonical basis on W n+l/ and divide the basis 

(/) • 1 J M • 1 

as e\ J = a, i = 1, . . . , n and e\ = e n+ij i = l,...,v. 
Now introduce on M. n+U a multiplication defined by 

„(/)>») _ Jm) Jf) _ i / V n „(/) (/) , v , p(m) (m) 

^ ' (/) (/) n ■ • i M M n 7 7 1 

e i e } = ° 3 i, J = !,-■-, n; el>e\>=Q, k,l = l,... 



Thus we identify the coefficients of bisexual inheritance as the structure 
constants of an algebra, i.e. a bilinear mapping of M. n+U x W l+U to M n+U . 

The general formula for the multiplication is the extension of (I3.ip by 
bilinearity, i.e. for z,t €G M. n+U , 

n v n v 
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using (I3.ip . we obtain 

*t = \ ELi (££=1 Ei=i p^C^ + «*»,•)) ei /} + 



(3.2) 



From (|3.2p and using (|2.2|) . in the particular case that z = t, i.e. x 
and y = v, we obtain 



« = ^ = ELi fe =1 EJ=i i^w) 4 /} + 



for any z G 5. 

This algebraic interpretation is very useful. For example, a BP state 
z = (x,y) is an equilibrium (fixed point, V{z) = z) precisely when z is an 
idempotent element of the set 5. 

If we write for the power (• • • (z 2 ) 2 • • • ) (t times) with z^ = z then 
the trajectory with initial state z is V t (z) = zM. 

The algebra B = By generated by the evolution operator V (see (|2.2|) ) is 
called the evolution algebra of the bisexual population (EABP). 

Remark 3.1. 1. If a population is free then the male and female types are 
identical and, in particular n = v, the inheritance coefficients are the same 
for male and female offsprings, i.e. 

p _ p(/) _ p( m ) 

The evolution algebra A associated with the free population is introduced 
and studied in Note that this algebra is commutative when the condition 
of symmetry Pikj = Puj *s satisfied, but it is not in general associative. In 
the next section we show that algebra B of bisexual population is commutative 
without any symmetry condition. Hence the algebra A is a particular case 
of the algebra B. 

2. It is easy to see that the EA introduced in is completely different 
from our EABP, i.e. of B. 

3. The algebra B is a natural generalization of a zygotic algebra for sex 
linked inheritance (see /5[ f77l[I^|[l^/ ). 



4. Basic properties of the EABP 
The following theorem gives basic properties of the EABP. 

Theorem 4.1. 1) Algebra B is not associative, in general. 

2) Algebra B is commutative, flexible. 

3) B is not power-associative, in general. 
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Proof. 1) Take e/\ ej such that P>J^ S ^ for some s and take such 
( f) 

that Pgi' 7^ for some r then 



(/) J m K c ( m ) 



i n 

-VP (/) e( / 'e W 



9=1 



7- P i/i- P ifc,t e r^' ) + non-negative terms 7^ 0. 



But 



J (/)/J" i )„H\ 

V e 7 e fc i 



0, i.e. fe^e^^e^/ePfe^ei" 1 ^ 



e i l c j c fc 7 - u ' \Pi J^k t- H 

2) Commutativity of B follows from formula (|3,2p . An algebra is called 
flexible if it satisfies z(tz) = (zt)z for any z,t. It is easy to see that a 
commutative algebra is flexible. 

3) To show that B is not a power-associative, in general, we shall construct 
an example of z such that (zz)(zz) ^ {{zz)z)z. Consider n = 1, v = 2. In 
this case 



11,1 



12,1 



1 pM , P ( m ) 
x i 11,1 ^11,2 



1 P H 1 P H 

x > ^12,1 ^12,2 



1. 



Denote a = p}7i> ^ = ^I2i - Take z = + e^ 71 . Then we have 

2 (/) 1 ("0 1 / -1 \ (rn) 

z = e\ + ae\ + [1 — a)e 2 ■ 



(4.1) z 2 z 2 



e[ f) + (a 2 + (l 



)b) e { ™ ] + (l-a)(a- 6 + 1)4' 



(m) 



Z 2 



+ J (a 2 + (1 - a)6 + a) + U\ - a) (a - b + 2)e 



(m) 



(m) 



(z 2 z)z = + i (a(a - b) 2 + (a + b) (a - b + 2)) 4 mj + 



(4.2) 



(m) 



|(1 -a) (3 + (a-6 + l) 2 )e 



(m) 



Assume a = P^l / 1 and a / 6 = P-^. Then from gj]) and flM]) we get 

(zz)(zz) {{zz)z)z. □ 



(m) 



5. B IS NOT A BARIC ALGEBRA 

A character for an algebra A is a nonzero multiplicative linear form on A, 
that is, a nonzero algebra homomorphism from A to R [13]. Not every alge- 
bra admits a character. For example, an algebra with the zero multiplication 
has no character. A pair (A, a) consisting of an algebra A and a character 
cr on A is called a baric algebra. In [13] for the EA of a free population it 
is proven that there is a character o~(x) = Yli x ii therefore that algebra is 
baric. But the following theorem says that the EABP, i.e. B is not baric. 
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Theorem 5.1. The EABP, B, has no a nonzero character. 

Proof. Assume a{z) = X^ILi a i x i + SJ=i bjVj, z = (x, y) G B is a character. 
We shall prove that a(z) = 0. For z = (x,y), t = (u, v) we have 

1 n v I n v \ 

"•(**) =JEE E %>42 + E WHS + 

i=l j=l \p=l q=l J 

n n n v v v 

a(z)a(t) = OiajXiUj + (i,b / (.r i i- / + it^-) + bibjyjVj. 

i=i j=i i=i j=i i=i j=i 

From cr(jrf) = cr(z)cx(t) we get 

ajOj = for any i, j = 1, . . . , n 

bibj = for any i,j = 1, . . . , v 

This system of equations has only the solution a\ = • • • = a n = b\ = ■ ■ ■ = 
b u = 0. □ 

Remark 5.2. In J13\j the baric EA of a free population is classified as al- 
gebra induced by a linear operator; unit algebra; constant algebra; Bernstein 
(stationary) algebra; genetic algebra; train algebra, etc. But since the EABP 
is not baric there are not similar algebras for bisexual population. 

6. B IS A DIBARIC ALGEBRA 

By Theorem 15.11 the algebra B is not a baric algebra. To overcome such 
complication, Etherington [5] for a zygotic algebra of sex linked inheritance 
introduced the idea of treating the male and female components of a popu- 
lation separately. In [11] Holgate formalized this concept by introducing sex 
differentiation algebras and a generalization of baric algebras called dibaric 
algebras. In this section we shall prove that the algebra 6 is a dibaric 
algebra. 

Definition 6.1. fl4\[2U^ Let 21 = (w,m)^ denote a two dimensional com- 
mutative algebra over R with multiplicative table 

9 9 / \ 

w = m = 0, wm = —{w + m) . 

Then 21 is called the sex differentiation algebra. 

As usual, a subalgebra B of an algebra A is a subspace which is closed 
under multiplication. A subspace B is an ideal if it is closed under multipli- 
cation by all elements in B. For example, the square of the algebra: 

A 2 = spanjzt : z,t E A} 

is an ideal. 

It is clear that 2t 2 = (w + m}^ is an ideal of 21 which is isomorphic to 
the field K. Hence the algebra 2l 2 is a baric algebra. Now we can define 
Holgate's generalization of a baric algebra. 
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Definition 6.2. An algebra is called dibaric if it admits a homomor- 
phism onto the sex differentiation algebra 21. 

Theorem 6.3. The algebra B is dibaric. 

Proof. Consider mapping <p : B — ¥ 21 defined by 

(6.1) <p(e"P)=w, i = l,...,n; ip(e^) = m, i = l,...,v. 

For z,t e R n+U , 



n v n v 

(x, y) = J2 xi( i f) + Yl %'4 m) ' t = ( u > v ) = Yl uie i f) + Yl v A 

i=l j=l i=l j=l 



using linearity of (p, f)6. 1[) . f)2. 1 1) and (j3.2|) we get 



n u /l \ 

1=1 7=1 ^ ' 



We also have 



(n i/ \ / n v 

Xi,w + ^2 yj m Yj Uiw + ^j 771 
i=l j=l ) \i=l 3=1 

i=l j=l V ' 

i.e. (p is a homomorphism. For arbitrary u = aw + f3m £ 21 it is easy to see 

that (p(z) = u if z = Yli=i x i e i + Sj=i 2/.7 e j e ^ with S=l ^ = a anc * 
X]J=i = Therefore (p is onto. □ 

Proposition 6.4. 17 i]/ an algebra A is dibaric, then A 2 is baric. 

As a corollary of this Proposition and Theorem 16.31 we have 

Corollary 6.5. The subalgebra B 2 is a baric algebra. 

Since B 2 is a baric algebra all types of particular algebras mentioned in 
Remark 15.21 can be defined and studied for the algebra B 2 . 

7. B IS A Banach algebra 
Define a norm || • || : B —> M. as follows 

n v 
i=l 3=1 

where z = Yh=i x i e \ f) + Y!j=i Vi e j G S - 

For a fixed a € £> consider the operator L a : B B, left (right) multipli- 
cation, defined as 

^a(z) = az (= za). 
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It is easy to see that the set < L (/) , L ( m ) , i = 1, . . . , n, j = 1, .., v > spans a 

I e * 6 J J 
linear space which is the set of all the operators of the left (right) multipli- 
cations. 

Proposition 7.1. Operator L a is a bounded linear operator for any a = 
Proof. For z = Ya=i x i e \ + Y^j=i yj e j £ B we have 

n n v \ 

= E h EE + ^) 4 S) + 

p=l \ i=l j=l ) 
9 =1 \ i=l j=l / 



n 

M*)ll = 5£ 



p=l 



i=i i=i 



+?E 



9=1 



i=l j=l 



< 



P=l \i=l j=l i=l i=l / 

^ v I n v n v \ 

9=1 \i=l j = l i=l j=l ) 

Using (I2.ip from the last inequality we get 

11^(^)11 < ||o||||y|| + W\\x\ < max{||a||, ||&||}||^||, 

where 



Ei a *i' \\ b \\ = J2\bj\, 

3=1 



X 



i=l 



Ew> imi = Ei%'I- 

i=l j=l 



Thus L a is bounded for any fixed a S B. 



□ 



Theorem 7.2. T/ie algebra B is complete as a normed linear space. That 
is, B is a Banach space. 

Proof. If {z k } converges then its all coordinates also converge and therefore 
limit of z k also will be an element of B. This completes the proof. □ 

Corollary 7.3. The algebra B is a non associative Banach algebra. 



10 



M. LADRA, U.A. ROZIKOV 



8. The derivations of B 

There are many papers devoted on the subject of derivations of genetic 
algebras (see e.g. [2,10,12,19]). In [12] an explanation of the genetic meaning 
of a derivation of a genetic algebra is given. For any algebra, it is known that 
the space of its derivations is a Lie algebra. The Lie algebra of derivations of 
a given algebra is an important tool for studying its structure, particularly 
in the non-associative case; so this is a natural development. 

In this section we describe the set of all derivations of B. Let D G Der(£?) 
and suppose 



D(e[ f) ) = E 4!^ + E 4^ » k = l,...,n. 

i=i i=i 



Die^) = E C'eP + E k = !,...,». 

i=l 1=1 

By the definition of derivation D(zt) = D(z)t + zD(t), we have 
D(4«ef) = D(e«)em+e«D(ef) = 



E 4I4 f) + E 4M" 4" + 4" E 4/^* + E &4 

\i=l 1=1 I \i=l 1=1 



1=1 \l=l 



^EfE(4 m ^ ) +4 m < ) ))^ ) = 

q=l \l=l J 

Consequently, 

(8.i) ; \ 

Z^=i ^"fci + "j; •'jw.g J - u ' 3) * i, . . . , 7i, g — i, . . . , v. 

Similarly, from Z)(ejj, ej ) = -D(e[ m ' ) )ej m ^ + D{e^) = we obtain 
£?=l (C^S + = °> s = 1, . . . ,n; j,k = 1, . . . ,v; 

(8.2) 



V™ (ff'fpW + ,/ m /pM\ _ n n 1c ft 



!,...,!/. 
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The equality D{e^ ) ef ii ) = D{e\ f) )e { ™ ] + ej /) D(eH) gives the following 
conditions 

Z^p=l ^"ip r pj,s a P s r ij,p ) 

ZU fe mp S - <s fp t<!) = °> M = 1, • • • , n; i = 1, . . . , i/; 

(8.3) 

( AH - rlf m F>^\ -L 

l^p=l y^ip r pj,t W pt r ij,p J "T 

Since for any z = (x,y) & B we have -D(-z) = J^™ =1 £i-D(e^)+5^ =1 Vj^i^^) 
D is uniquely defined by the matrix V = V{D) = (d{J , d{™, d™J , d. 
Hence, 



mm 

ij T~vp '"'pi '"'pq / y=i n 

p,q — 1 , .. .,W 



Der(£) = {£> : V(D) satisfies JS3J, ([821), §8^} ■ 

(f) (f) 

Example. Consider the case n = 1, v = 2. In this case Pq ^ = = 1. 
Denote pf™} = a, pf^ = 1 - a, P-f™} = 6, pf^ = 1-6- Assume a = 1 - b, 
a^l. Then using ([OP . we get 

£)( e C«) = 0; D(e { ™ ] ) = a ■ (e< m) - e { fy; D(e { ™ ] ) = a • (-e£ m) + ej, m) ), 

where a G R. Consequently, for arbitrary element z = (£,2/1,2/2) € R 1+2 we 
have 

Z^) = a.( yi -y 2 )(4 m) -4 m) ). 
Thus if n = 1, z/ = 2 and a = 1 — b, then 

Der(£) = [D : D(z) = a( Vl - y 2 )(e[ m) - e^), z = (x, yi ,y 2 ) G B,a G r} . 

9. Dynamics of the operator ( 12.21) 
Extend the operator (|2.2|) on M. n+U , i.e. consider the operator V : R n+I/ — > 



z = (x,y) z' = (x',y') = V{z) defined as 

(9-i) = E 45^ ^' = E p SW 

i,fe=l i,fe=l 

Consider the following linear form X : R n — > R ( Y : W — >• R) defined by 

n / u \ 

(9-2) X(x)=E^' k(y) = E^ " 

i=i V fc=i / 

Denote 

(9.3) ^ = { z =( x ,y):X(x)=y( 2 /) = i}, i = 0, 1, 

the product of the z-hyperplanes in W 1 and M 1 ", respectively. 
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A point z = (x, y) € IR n+iy is called a fixed point (resp. zero point) of V if 
V(z) = z (resp. V{z) = 0). 

Proposition 9.1. (1) If z is a fixed point then z G i^o U 
(2) If z is a zero point then z E {z = (x,y) : X(x)Y(y) = 0}. 



Proof. From (J9JJ) and using (J2JJ) we get 

x(x') = E7=i *} = £5=i (£"=i 43) *w* = 

= £r=i = £IJ =1 (£r=i 4?) *w* = x{x)Y(y). 

(1) If z is a fixed point then X(x') = X{x) and Y(y') = Y(y), this by 
(19~3D gives that X(x) = Y(y) and X(x) = (X(x)) 2 . Hence X(x) = 
Y(y) = or 1. 

(2) If z is a zero point then X(x') = Y(y') = 0, and gives X(x)Y(y) = 
0. 

□ 

Let z^ = (x(°),?/ )) E M n+l/ be an initial point. Its trajectory is defined 
by *(*) = (arW,yW) = U^*" 1 )), t = 1,2,.... Denote by the set of 

limit points of the trajectory {z^}^ . If uj(z^) consists of a single point, 
then the trajectory converges, and oj(z^) is a fixed point of the operator 
V. 

The following theorem gives an upper estimate of the set oj(z^) 



Theorem 9.2. We have 



w(z (0) ) C { 



' H if \X(xW)Y(y(°))\ < 1 



H x if \X{x^)Y(y(°y)\ 



1 



where H n 



if |X(x(°))r(y(°))|>l, 
{z = (x,y) : X{x) = Y{y) = +oo}. 



Proof. Using (@3J we get X(x^) = Y{y®) = X{x^)Y{y^) for any 
t = 1,2,.... Iterating this recurrent equation we obtain 

(9.5) X(x®) = Y(y®) = (X(x(°))y(y(°))) 2t * . 

Consequently, 

' if \X(x<®)Y(y<®)\ < 1 

lim X(x {t) ) = lim Y(y W ) = < 



t— >oo 



t— >oo 



1 if |X(x(°))Y(y(°))| = 1 
+oo if \X{x^)Y(y^)\ > 1. 



□ 
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10. A SPECIAL CASE OF AN EABP 

In this section we consider a special case of an EABP giving an addi- 
tional condition on heredity coefficients (|2.1j) . i.e. consider the coefficients 
as follows 



(10.1) P 



(/) 

ik,j 



'''.1 



if k = 1 



1 if k^l,j = l 

I if k^lj^l. 



P 



(m) 



ikl 



( b k i if i = 1 

1 if i l f l = 1 

{ if 1,1^ 1. 



The matrices A = (aij) and -B = (bki) by (12. ip satisfy the following condi- 
tions 
(10.2) 

n v 

aij>0, y^Qjj = 1, i = l, ...,n; fefc/ > 0, y^bfc; = 1, k = l,...,v, 
j=i i=i 

i.e. both matrices are stochastic. 

Remark 10.1. The condition $10. 1\) is taken to simplify our computations; 
in this way we reduced both cubic matrices ^-Pf/^Jj (Pikfj t° quadratic 
matrices. But that condition has very clear biological treatment: the type 
"1" of females and the type "1" of males have preference, i.e. any type of 
female (male) can be born if its father (mother) has type "1". If the father 
(mother) has type ^ 1 then only type "1" female (male) can be born. 



Under condition (jlO.ip the multiplication (|3.ip became as 
(10.3) 



( EU^+ZUhie™, if i = l,fc = l 



e i e k 



Sf) 



[ m) , if i^l,k = l 



e\ J ' +Y,Uhie\ m> , if i = l,k^l 



e { pe { p = Q, i,j 



I.. 



e[ /} +e { ™\ if i + l,k + 1. 

(m) (m) 

,n; e y k > e\ 



0, k,l 



1, 



Operator ()2.2[) has the following form 



(10.4) 



E"=i (a.12/1 + EL2 2/*)^ 

yi E"=i a 'i x i' j ' ^ 1 

ELi ( 6 fcisi + E™=2 z/fc 
y[ = x iJ2k=i b kiyk, i^i 

Denote by B\ the EABP defined by the multiplication table f j 1 . 3 1) . 
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10.2. Idempotent elements of B\. A element z G B is called idempotent 
if z 2 = z; such points of an EABP are especially important, because they 
are the fixed points of the evolution map V, i.e. V{z) = z. We denote by 
Id(B) the idempotent elements of an algebra B. Clearly, £ Id(B) and this 
set is an algebraic variety. By Proposition 19.11 we have Id(B) C Hq U H\. 
In this subsection we shall describe idempotent elements of B\. First we 
describe the set of idempotent elements which belong in Hq Dld(Bi) and 
after that we shall describe elements of H\ DZd(Bi). 
Idempotents in Hq. 

Using (|10.4p and the condition that z € Hq, from the equation V{z) = 
z 2 = z we obtain 

xj = yi E^Oij - aij)xi, j = 2, . . . , n; 

(10. 5) 

yi = x iYlk=2( b ki - hi)yk, i = 2,...,u. 

Case xiyi = 0: If x\ = then yi = ■ ■ ■ = y v = 0, consequently, y\ = 0. 
Similarly, if y\ = we get x\ = ■ ■ ■ = x n = 0. Hence in case xiyi = we 
have unique idempotent z = 0. 

Casexiyi ^ 0: Consider matrices C y = (Qj)jj = 2....,n and D x = (dki)k,i=2,..., 
such that 

^ J (ojj - a Xj )y if i / j 

[ (Oy - aijjy - 1 if z = j 

Then equation (|10.5p can be written as 

(10.6) C yi x = 0, D xi y = 0, 

where x = (x 2 , . . . ,x n ), y = (y 2 , ■ ■ .,y v ). 

Consider now x\ & M. \ {0} as a parameter, then equation D xi y = has 
a unique solution y = if det(D Xl ) ^ 0, which gives y\ = 0, i.e. this is a 
contradiction to the assumption that y\ ^ 0. 

If det(D Xl ) = then we fix a solution x\ = x\ ^ of the equation 
det(D Xl ) = 0. In this case there are infinitely many solutions y* = (y|, . . . , y*) 
of -D^y = 0. Substituting the solution y* = — Ylk=2 ul m ^yi x = ^, we get 

n 

(10.7) = yl y^( a ij ~ a ij) x i = 0, j = 2, . . . ,n. 

i=2 

This system has a unique solution xi = ■ • • = x n = if det(C y *) 7^ but in 
this case we get xi = which is a contradiction to the assumption that x\ 7^ 
0. If det(C y *) = then we have infinitely many solutions x* = (x^, 
Hence we have proved the following 

Proposition 10.3. We have 

HQC\Zd(Bi) = {0}U 

{((xl,...,x* n ),(yl,...,y:)) : C v *x* = 0, D x *y* = 0,det(L»^) = det(C y j) = 0}. 
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Idempotents in H\. Using (I10.4P and the condition that z G Hi, the 
equation V(z) = z 2 = z can be written as 

xi = YJt =l {{an-l)yi + l)xi\ 
x j = 2/1 2^=1 a ij X i' 3 T 1 - 

Consider several cases. 

Case x\=yx = 0. In this case we get z = which is not in H±. 

Case x\ 7^ 0, y% = 0. In this case xi = 1, xj = 0, j = 2, ...,n. 
Substituting xi = 1 in the system of equations f)10.9j) we get 

(io.io) = s ; 2te 

yi = Hk=2 b kiyk, l = 2,...,u. 

Denote £ 2 = {b' kl ) kl=2 ^ with 

, J bki-1 \i k = l 
{ b M if k ^ I. 

If det(i? 2 ) 7^ then B 2 y = gives y\ = ■ ■ ■ = y v = but this does not sat- 
isfy Ylk=2 Vk = 1- ^ det(Z? 2 ) = then = has infinitely many solutions 
y = (y 2 , • • • j )• We must take these solutions which satisfy Ylk=2 bkiVk = 0. 
So in this case we have the following idempotent elements (which belong to 
Hi) 

_ f {((l,0,...,0),(0,y 2 ,...,^)):i?2y = 0,EL 2 ^iyfc = 0} if det(S 2 ) = 
° ~ \ if det(-B 2 ) ^ 0. 

Case xi = 0, y\ 7^ 0. This case is similar to the previous case. Denote 
A 2 = ( a'„- ) , with 



jj'=2,...,n 

, / - 1 if « = j 
Ojj if i / j. 



h 



Then we have the following idempotent elements (which belong to H\): 

f {((0, x x , . . . , x n ), (1, 0, . . . , 0)) : A 2 x = 0, £™ =2 c^x* = 0} if det(A 2 ) = 

~\ if det(A 2 ) / 0. 

Case xi 7^ 0, y\ 7^ 0. Here, to avoid many special cases we assume 
that an 7^ 1- Take y\ as a parameter and solve the system (|10.8p which is 
equivalent to the following 

(10.11) E( fl ^i: fl ^ 1+1) +^)^=^ ^ n. 
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Denote U y = {uij)ij=2,...,n with 



aij(Ki-i)y+i) 



+ aijy if i^j 



c kl 



1-ou ^ y 
" 1 + «, y y - 1 if i = j 

Subcase det(U y ) ^ 0. In this case we have x% = ■ ■ ■ = x n = 0, conse- 
quently, x\ = 1. For x\ = \ from the equation (|1U.9|) we get = where 
B i = ( b ll)k,l=i,...,u, with 

( b k i -1 il k = l 

b h = s 

[ b H if k^l. 

So in this case we have the following set of idempotent elements of Hi. 

_ f {((l,0,...,0),(yi,...,y l/ )):Siy = 0,det(t^)^0}if det(Bi) = 

2 ( « if det(Si) + 0. 

Subcase det(C/ 2/ ) = 0. In this case we fix a solution yi = y* of det(U y ) = 
0. We have infinitely many solutions (x^, Denote C x = (c* kl )k,i=i,...,v, 

with 

' (&n - l)x if k = 1 = 1; 
(6 fe i-l)x + l if = 1; 

fy/x — 1 if A; = I 7^ 1; 
ftfcix if k ^ 1,1 = 2, . . . 

Now one has the following set of idempotent elements 

h = {((x*i,...,x* n ),(y*i,...,yt)):C xt y* = 0, U yt x* = 0, 

det(C^) = det(O 8 j) = 0}. 
Thus we have the following 

Theorem 10.4. // an / 1 £/ien £/ie /u/Z sei of the idempotent elements 
which belong to Hi is 

Xd(Bi) n Hi = I U h U I 2 U J 3 . 

Remark 10.5. 5y Proposition \9.1\ we can conclude that Proposition \10.^ 
and Theorem TJTj^give the full set Xd{Bi). But conditions described in the 
Proposition and Theorem are complicated in general. In the sequel of this 
subsection we shall consider additional conditions on {flij), (b«) and under 
these conditions we explicitly describe the setld(Bi). 

Now we consider a particular case and describe the full set of idempotent 
elements of £>i. 

If we consider the case 

{1, % = j \ 1, k = I 

bki 
0, i + j, 

then the following is true 
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Proposition 10.6. If the condition \1Q.12§ is satisfied then 
Xd(Bi) = {0}uU(l,x 2 ,...,x n ),(l,y 2 ,.. .,?/„)) : = ^y k = olu 

I i=2 k=2 J 

< ((l,x 2 , . . .,x n ), (1,3/2, • • - ,2/1/)) : = = _1 f U 

|((l,0,...,0),(yi,...,^)) :^y Jfe = l,y 1 ^l|u 

. . . ,x n ), (1,0, . ..,0)) : = l,xi / l| . 

Proof. The equation = z, for z = (x,y) £ M n+I/ , using condition ()10. 12[) 
can be written as 

(1 - yijxi = Yl=2 Vk E?=i x i> 
(l-yi)zj=0, j = 2, ...,n 

(10.13) 

(i - = Ya=2 x i ELi 

(l-xi)yz = 0, Z = 2, . 
The simple analysis of the system ()10.13p gives the set of all idempotents. □ 

10.7. Absolute nilpotent elements of B\. The element z is called an 
absolute nilpotent if z 2 = 0, i.e. z is zero-point of the evolution operator V. 
For £>i the equation z 2 = is equivalent to the following system of quadratic 
equations 

E"=i Kiyi + EL2 J/*) x * = °; 

2/1 Ya=i = °> i = 2, • • • , n; 
(10.14) _ 

x i Y!k=i hlVk = 0, / = 2, . . . , v. 

By Proposition 19. II we have ElLi x « Efc=i 2/fc = 0- There are the following 
three cases. 

• Case EILi x i = Ylk=i Vk = 0- I n this case, from the system of 
quadratic equations (jlO. 14|) we get 

Vi J2i=i a *i x i = °j i = !> ■ ■ ■ > n ; 

(10.15) 

x i Efe=i b fci2/fe = 0, 1 = 1,..., v. 
Subcase y\ = x\ = 0. In this case one easily gets the following set 
of solutions 

-^00 = \ ((°) x 2, ■ ■ ■ , x n), (0, 2/2, ■ ■ ■ , J/j/)) : ^ »i = 0, ^2 Vk = \ ■ 

[ i=2 k=2 J 
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Subcase x\ = 0, y\ ^ 0. Denote Aq = (ajj)i,i=2,...,n- It is easy to 
see that the set of solutions is as follows 



- / vm 



'01 

where 



A" 1 if det(Ao) ± 0; 
A" if det(Ao) = 0, 



AA 1 = l((0,...,0),( yi ,...,y u )) : yi ^0, j^y fc = o|, 

f n n 

A/" = < ((0,x 2 ,. • • ,x n ), (yi,. . . ,y„)) : ^x^ = ^ a a x; = 0, 



i=2 i=2 



A)X = 0, j/i / 0, y fc = ^ . 

k=l J 

Subcase x\ ^ 0,y\ = 0. Denote i?o = (pkl)k,l=2,...,v In this case 
the set of solutions is 

' M 2 if det(B ) f 0; 
A/" 3 if det(5 ) = 0, 



A^ 



'10 

where 

( n •> 



i=l 

n 



N 3 = < ((xi, . . . ,x n ), (0,1/2, • • • ,y v )) : xi / 0,^Xj = 0, 
I i=i 

Boy = 0, ^ yk = ^2 bklVk = \ ■ 

k=2 k=2 ) 

Subcase x\ ^ 0, y\ ^ 0. In this case it is easy to get the following 



if det(^) / or det(S) + 0; 
A/" 4 if det(A) = det(B) = 0, 



<i = 
where 

M 4 = {((xi, . . . , x n ), (yi, . . . , y,)) : Ac = 0, By = 0} . 

• Case x « = 0' Sfc=i 7^ 0- I n this case the system of quadratic 
equations (|10.14p can be written as 

YZ=i (^i^ 1 + T!i=2 Vk) xi = o 

Vi Ya=i a ij x i = °, j = 2,...,n 

(10.16) 

x i ELi - !) Vk = o 

x i Z)fc=i hiVk = 0, l = 2,...,u. 
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Subcase y± = x% = 0. In this case we have the following set of 
absolute nilpotents 



•A/qo = \ ((0. x 2, ■ ■ ■ , x n ), (0, y 2 , • • • , 2/iz)) : x < = °' Yl Vk ^ 

I i=2 k=2 

Subcase x\ = 0, y\ ^ 0. From (jlQ. 16f> we get 

E"=2 i a HVl + Efc=2 2/fc) ^ = o 
( 10 - 17 ) ElLi^i = 0, j = 2,...,n. 



This has the following set of solutions 



A" 5 if det(Ao) / 0; 
A" 6 if det(Ao) = 0, 



where 

AA 5 = ( ((0, . . . ,0), (2/1, . . . ,y u )) : 2/1 # 0, ^yfe + 



fc=i 

n n / 



A/" 6 = < ((0,x 2 ,.. • ,x n ), (2/1,. • .,2/1/)) : J2 Xi = Y [ ailVl + ^2 Vk ) Xi = °' 

[ i=2 i=2 \ fe=2 / 

i4ox = 0,yi^0,^y fc ^ol. 

fc=i J 

Subcase x% 7^ 0, 2/1 = 0. In this case from ()10.16p we get 

(10.18) EiU(o«-D»=o 

Efc=2 & fc«2/fc = °> l = 2,...,u. 

If det(So) / we get 2/2 = • • • = Vv = 0, this is impossible since 
we have condition EaUi 2/fc 7^ 0- For det(-Bo) = the set of absolute 
nilpotents will be 



■A/10 = S ((^l)- • • > x «)' ( ,2/2,- ■ ■ ,2/1/)) : xi ± 0, = 0, 

I i=i 

B y = 0,J2 Vk + 0, ( 6 « - 1) 2/fc = L 

fc=2 fc=2 J 

Subcase x\ ^ 0,y\ ^ 0. Denote A = (ay)^-^...^, with 
J 1 if i = 1; 
1 Oij if j ^ 1. 
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and B = (b k i)k,l=i,...,u, with 

In this case we have 
Mil = 



bfci - 1 if / = 1 
b k i if I + 1. 



j if det(A) / or det(B) ^ 0; 
[ AT 9 if det(A) = det(B) = 0, 
where 

TV 9 = {{(xi, x n ), {yi, . . . , y v )) : Ax = 0, By = 0, 

n / ^ \ v \ 

^2 [ aiiyi + yfc ) Xi = °' X] yfc ^ r • 

i=l V fc=2 / fc=l J 

• Case Y^i=i Xi ^ ^' Sfe=i ^fc = 0- This case is similar to the previous 
case, to get the set of nilpotents, one has to change sets A/^ , i, j = 0, 1 
as follows. Replace x and y, also rearrange parameters with 
Denote resulting sets by Mfj, i,j = 0, 1 respectively. 

Thus we have proved the following theorem 

Theorem 10.8. The full set TV of absolute nil-potent elements of the algebra 
B\ is 

2 
i=0 
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